Integral Representations for Free Energies of Macroionic Suspensions and Equation of 
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A generating functional which results in the Poisson-Boltzmann equation and boundary condi- 
tions for an average electric potential of a macroionic suspension through an extremal condition 
is constructed in a mean field theory. The extremum of the generating functional turns out to be 
identical with the Helmholtz free energy of the system which has an integral representation in terms 
of the average electric potential satisfying the Poisson-Boltzmann equation. From the Helmholtz 
free energy, the chemical potentials of small ions and chemical potentials of effective valencies of 
macroions are calculated and, as the total sum of them, an integral representation of the Gibbs free 
energy of the system is derived. Difference of two free energies leads to an equation of state for 
osmotic pressure of small ion gas in an environment of macroions in the suspension. 
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In the Debye-Hiickel theory [1] which is a prototype 
of mean field description, the Helmholtz and Gibbs free 
energies of an electrolyte, F and G, are calculated in a 
self-consistent manner by linearizing the PB (Poisson- 
Boltzmann) equation. The difference G—F being identi- 
cal with the equation of state for osmotic pressure turns 
out to consist of a term of the ideal gas law, i.e., the van't 
HofF law, and a negative-definite term representing the ef- 
fect of electric charges. Interpreting this decrease of the 
osmotic pressure as a manifestation of attractive nature 
of a net electric interaction and applying the DH (Debye- 
Hiickel) theory to macroionic systems, Langmuir [2] rec- 
ognized first a possible existence of a long-range Coulom- 
bic attraction between macromolecules with like electric 
charges and predicted that the equation of state resem- 
bling a van der Waals isotherm leads to an instability and 
phase separation of ionic species. Verwey and Overbeek 
had criticized strongly this ingenious insight by pointing 
out "the danger involved in the use of the linear approx- 
imation characteristic of the DH theory" [3] . 

Recent development of experimental techniques has re- 
vealed existence of inhomogeneities of particle distribu- 
tion in colloidal suspensions such as clusters and voids [4- 
7]. These new large scale phenomena, along with well- 
known observations as salt-induced-melting of colloidal 
crystals [8], crystalization through multi-stage phase 
transitions [9] and reentrant solid-liquid transitions [10], 
stimulate strongly studies of charge-stabilized macroionic 
systems. Therefore, it is natural and reasonable to revisit 
the long-ignored idea of Langmuir which had predicted 
an occurence of phase separation in macroionic suspen- 
sions. The purpose of this article is to formulate a mean 
field theory for charge-stabilized suspensions of macro- 
molecules [11, 12] to revive the Langmuir's viewpoint 
without recourse to a linear approximation. 

Direct application of the results of the DH theory to 
macroionic systems in the original Langmuir article [2] 
led necessarily to an unphysical ionic democracy [12] 



where all the ionic species are treated so symmetri- 
cally that small ions must have surrounding clouds of 
macroions as well as macroions have clouds of small ions 
and small ions have mutual atmospheres of small ions 
with opposite charges. To remedy this apparent draw- 
back of Langmuir's treatment, we adopt a viewpoint of 
ionic aristocracy [12] where the size and the time scale 
of motion of macroions are considered to be asymmet- 
rically larger than those of small ions. Postulating that 
macroions are adiabatically fixed, we derive integral rep- 
resentations of the Helmholtz and Gibbs free energies of 
macroionic suspension F and G in terms of an average 
electric potential satisfying the PB equation. The dif- 
ference G — F leads to the equation of state for the os- 
motic pressure of small ion gas in a fixed configuration 
of macroions. 

We investigate a macroionic suspension with M num- 
bers of macromolecules in a container with a sufficiently 
large volume fi. Henceforth macromolecules are called 
particles for brevity. Effective values of the surface 
charge [10], volume and surface area of the n-th parti- 
cle are Zne, w„ and Sn, respectively. The gas of small 
ions with valency Zj occupy the volume 



V 



M 



(1) 



As a basic postulate of the present scheme, the thermody- 
namic function Q of the system, viz., the internal energy 
U, the Helmholtz free energy F and the Gibbs free en- 
ergy G, is generically postulated to consist of the sum of 
its electric part and non-electric part Q'^ as 



(2) 



To make arguments clear, we model macroionic sus- 
pensions as follows: 

1. Solvent of the system forms a uniform medium of 
background with a dielectric constant e . 
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2. Particles take a static distribution. The config- 
uration of tlie particles is represented by the set 
{R} = {Ri,R2, - ■ ■ Rm} of their center-of-mass 
coordinates. 

3. In a static environment of particlcsCthc gas of small 
ions reaches rapidly to a thermal equilibrium in 
which the gas of j-th species of small ions with va- 
lency Zj {j — 1, - ■ ■ ,s) has a number density given 
by the Boltzmann distribution 

r^,■W = %„e-^^'^(^) (3) 

where 4>{r) is a ratio of the ordinary average elecrtic 
potential ip{r) to an average thermal energy kT = 
1/(3, i.e., 

<j){r) = f3eij{r). (4) 

The coefficient rij^ in Eq.(3) is a constant deter- 
mined by the number Nj of small ions in the vol- 
ume V as 

Nj=n,, [ e-^^'f^^UV. (5) 
Jv 

4. The dimensionless potential ^{r) satisfy the PB 
equation 

V^0(r) = -47rAB ^^,n,(r) (6) 

and the boundary conditions 

i^n ■ V(/)(r) = -47rAi3Z„f7„(r) (7) 

on particle surfaces Sn{n = 1, • • ■ , M) where A_b is 
the Bjerrum length Xb = e^fi/eGun is a unit vec- 
tor normal to the surface Sn and (T„(r) is a surface 
charge density normalized by the condition of sur- 
face integral 

/ a„(r)d5„ = l. (8) 

JSn 

For brevity, Eq.(7) is considered for n = to rep- 
resent the boundary condition of the inner surface 
of the container 5*0 which has a surface charge Zqc, 
a surface charge density CTo(r) and a unit normal 
vector i/Q. 

In analytical mechanics, the basic equations for vari- 
ous systems are described in a unified way by the Euler- 
Lagrange equation which is derived from the principle of 
minimum action. The macroionic suspensions are quite 
different from the mechanical systems. Nevertheless, it 
turns out possible to construct a generating fuctional 
from which the basic equation and the bounary condi- 
tions, viz., Eqs. (6) and (7), are deduced by imposing 



the extremal condition upon it. In factCwe can prove 
the following theorem. 

[Theorem 1] The generating functional for the PB 
equation and its boundary conditions is given by 

mx] = j^{^x{r)fdV 

M „ 

+ ^Zn(j) x{r)(Tn{r)dSn 

n=l •^■5" 

-^iV,ln(^l J^e-^^^^^UV^ (9) 

where x('") is an unknown smooth functionCthe volume 
integral is taken over V and the surface integrals are done 
on the surfaces of particles and the inner surface of the 
container. 

To examine the extremal condition of the functional, 
let us calculate the difference between ^[x-|-(5x] and J^[x] ■ 
Applying the Gauss integral theorem, we obtain 

5l3r = l3:F[x + 5x]-mx] 

+ -t\- V / (l^n • VX + ^■KXBZn<Jn)5xdSn. (10) 

For the functional ^[x] to take an extremum at xii") = 

(j){r)C5(3J-' must be an infinitesimal of second order with 
respect to an arbitrary variation dx as 

spj^ = m + ^x] - m<t>] = o5{x^). (11) 

This implies that the right hand side of Eq.(lO) which 
consists of infinitesimals of the first order must vanish at 
x(t) = 0(t). Therefore, the extremal condition of the 
functional T[x\ leads simultaneously to the PB equation 
(6) and the boundary conditions (7) by using the rela- 
tion between Eq.(3) and Eq.(5). For later arguments, 
it is convenient to express the extremal condition of the 
generating functional at x = </> in Eq.(ll) as 

5(3r[(l)]/5(l) = 0. (12) 

As shown below, the generating functional T is deeply 
related with the Hclmholtz free energy of the system. 

The electric part of the internal energy J/"' is identical 
to the average value of electric energy E of the system 
which is given, in the mean field description, by 

PU^' = pE=^^ j^{SI4>fdV 

= ^Z] / ZjrijHV+^^Znf (Aa„dS„.(13) 
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The right hand side of this expression is obtained by using 
the Gauss integral theorem and the PB equation with 
associated boundary conditions. This estimation for E is 
exact in the mean field description. 

The Helmholtz free energy F and the internal energy U 
of the system are related by the themodynamic relation 



V 9(3 J v,Nj,M,{R} 



(14) 



where the derivative with respect to f3 is taken for the 
fixed values of the volumeCthe numbers of small ions and 
particles and also for the fixed configuration of particles. 
In the present scheme, the same relations are assumed 
to hold separately both for the electric parts f ^ and If^^ 
and the non-electric parts and of the Helmholtz 
free energy and the internal energy. 

Here we use the thermodynamic relation in Eq.(14) for 
F^^ and If^^ = E to find the following theorem. 

[Theorem 2] The electric part of the Helmholtz free 
energy of the system F'^^ is the extremum of the gener- 
ating functional, i.e., F^^ = J^[4']. 

To prove this theorem we notice that J-[(t)\ depends on 
P through the Bjerrum length As and the average electric 
potential (p. It is readily shown that 

fdp^\ ^ d\B ( di3:F[ct>] \ dcj) ( 5(3:f[4>] \ 

\ d(3 ) d(3\ OXb J^'^dpy Sep 



= E 



(15) 



where the condition of extremum in Eq.(12) is used. As 
is readily confirmed it is possible to prove this theorem 
also by using the Debye charging-up formula. 

In this way the Helmholtz free energy of the macroionic 
suspension F is expressed in terms of the average electric 
potential (f) satisfying the PB equation as follows: 

F = + F° 

= -U'TJ^Zj f njcl>dV+hTj2Zn(f <t>CTndSn 

-kTj2Nj\n(^^ J e-'^^dV^ +F°. (16) 

Notice that this integral representaion is exact in the 
mean field description. If we use the identical relation 
in Eq.(13). the Helmholtz free energy F can take various 
different but equivalent forms being convenient for each 
purpose. 

Integration of the PB equation in Eq.(6) over the vol- 
ume V and usage of the Gauss integral theorem for the 
boundary condition in Eq.(7) result in the condition of 
electric neutrality 



(17) 



In the macroionic suspension where the number and con- 
figuration of particles are fixed, this condition requires 
that the effective valency of particles Z„ must change 
against a variation of the number of small ions Ni . There- 
fore the chemical potential of the small ions which is 
obtained by dF/dNj must necessarily be correlated to 
the derivative of F with respect to Z„. It is reasonable 
to interpret the derivative of the Helmholtz free energy 
dF/ dZn as the chemical potential of effective valency of 
the n-th particle [11, 12]. The Gibbs free energy of the 
system G is obatined by summing up the chemical po- 
tentials of the small ions of all kinds and the chemical 
potentials of the effective valencies of all particles. 

To take derivatives of F with respect to Nj and Z„ 
and to sum up the component chemical potentials with- 
out contradiction to the constraint condition [13], we in- 
troduce the following differential operator as 



D 



d 



NZ 



d 

dz: 



(18) 



It is readily confirmed that, when applied to the neu- 
trality condition, this operator reproduces the condition 
again without resulting in any additional constraint [13]. 
Therefore, this operator Dnz can be interpreted as a 
kind of covariant derivative which preserves invariance of 
a surface spanned by the constraint in a fictitious space 
(iVi,---,iV„Zi,---,ZM). 

For application of the covariant derivative D^z, it is 
convenient to use the Helmholtz free energy in the form 

PF = --^ [ {V</.(r)}2dF + VZ„ / (^{r)an{r)dSn 
SttAb Jv „ Js„ 

- E 1^ {y e-^^*(^)rfy^ + PFo. (19) 

Noting that the average potential as the solution of the 
PB equation depends on Nj and Zn, we find 



+Dnz{(P) 



5f3F \ 
H ) 



(20) 



Ni,Z„ 



Then, using the extremal condition in Eq.(12) and the 
definition G = Dnz{F), we obtain the theorem. 

[Theorem 3] The Gibbs free energy of the macroionic 
suspension has the integral representation 



G = kTY^Zuf 0(r)a„(r)d5„ 

-kT^Nj\n(^ j e-^'^^^UV^ + Go (21) 

where Gq = Dnz{Fo) is a non-electric part of the Gibbs 
free energy. 
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Among the three integral representaions in Eqs.(13), 
(19) and (21), tliere liolds the relation 

G-Gq=^F-Fq + E. (22) 

On the other hand, the difference between G and F is 
related to the osmotic pressure of the suspension as 

G - F = (Po + P)V (23) 

which reduces to 

Go - Fo = PoV + kTY, (24) 

3 

at the limit where all the solutes become neutral. 
HereCPo^ and kTJ2j^j are contributions from the 
pressure term of solvent and the van't Hoff term. There- 
fore, the gass of small ions in the environment of 
macroions is subjected to the equation of state 

PV = kTY,Nj +E (25) 

j 

which is exact in the mean field description. 

In this way we have derived the equation of state for 
osmotic pressure of small ions without resort to the lin- 
ear approximation. The right hand side of Eq.(25) in- 
cludes the average electric energy E{{R}). In contrast 
to the equation of state in the DH theory of a simple elec- 
trolyte where the net contribution of electric charges is 
a negative-definite constant, the average electric energy 
E{{R}) takes various values depending on the configu- 
rations of particles. Which sign does E take? It is safely 
argued that E is negative for almost all thermodynami- 
cally equilibrium states. To testify to this assertion, it is 
sufficient to apply the condition of electric neutrality to 
an arbitrary region with medium size in the suspension. 
Ions with unlike charges tend to come closer and ions 
with like charges tend to push farther apart in average 
so that the condition of electric neutrality is maintained 
in the region. As a result the net electric energy becomes 
negative. The electric energy may be instantaneously 
positive in a region in a suspension where particles take 
an extremely singular configuration. Such a configura- 
tion, however, is not thermodynamically stable. 

Therefore it is possible to conclude that, although the 
mathematical formulation of his theory was immature, 
the insight of Langmuir was correct. The simple form 
of Eq.(25) is a direct manifestation of the Langmuir's 
hypothesis of the Coulombic attraction. An electrically 
stable configuration of particles with like charges which is 
established through the intermediary of small ions among 
particles tends inevitably to attract and cage the small 
ions, and consequently reduces the osmotic pressure. The 
simple form of Eq.(25) enables us to proceed computer 
simulations. Notice that, if the equation of state includes, 
not the electric energy E, but a free energy F or G, 
its physical implication becomes much more indirect and 
computer simulations becomes difficult. 



In the mean field theory with the linearized PB equa- 
tion [11, 12], the equation of state in Eq.(25) had been de- 
rived and the approximate Gibbs free energies had been 
proved to have a medium-range strong repulsive part and 
a long-range attraction tail. Thermodynamic processes 
of particles develop slowly toward a direction in which 
the value of the Gibbs free energy G{{R}) decreases. It 
is the long-range attractive tail of the pair effective inter- 
action included in G{{R }) that is interpreted to promote 
the formation of large scale phenomena such as clusters 
and voids. In this connection it is worthwhile to point 
out that the Langmuir theory is related with volume- 
term theories [14, 15] which is now advocated to explain 
phase separations in macroionic suspensions by a spin- 
odal instability sensitive only by small ions. Schmitz et 
al. [16, 17] clarified mutual relationships among the Lang- 
muir theory, the theory with the Gibbs attractive inter- 
action and the volume-term theory essentially within the 
extent of linearized mean field description. The present 
formalism of mean field theory without linear approxi- 
mation strengthens validity of their arguments. 
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